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Abstract 



Abstract. In particle therapy, knowledge of the stopping-power ratio (STPR) of the 
ion beam for water and air is necessary for accurate ionization chamber dosimetry. 
Earlier work has investigated the STPR for pristine carbon ion beams, but here we 
expand the calculations to a range of ions (1 < z < 18) as well as spread out Bragg 
peaks (SOBPs) and provide a theoretical in-depth study with a special focus on the 
parameter regime relevant for particle therapy. 

The Monte Carlo transport code SHIELD-HIT is used to calculate complete 
particle-fluence spectra which are required for determining STPR according to the 
recommendations of the International Atomic Energy Agency (IAEA). 

The STPR at a depth d depends primarily on the average energy of the primary 
ions at d rather than on their charge z or absolute position in the medium. However, 
STPRs for different sets of stopping-power data for water and air recommended by the 
International Commission on Radiation Units & Measurements (ICRU) are compared, 
including also the recently revised data for water, yielding deviations up to 2% in the 
plateau region. In comparison, the influence of the secondary particle spectra on the 
STPR is about two orders of magnitude smaller in the whole region up till the practical 
range. The gained insights enable us to propose simple analytical expressions for the 
STPR for both pristine and SOBPs as a function of penetration depth depending 
paramctrically on the practical range. 



PACS numbers: 87.55. Qr, 34.50.Bw, 87.53.Bn, 87.55.K- 
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1. Introduction 

Stopping powers are essential for calculating the dose deposited by ionizing particles. 
The deposited dose is described as the mass stopping power multiplied with the 
particle fluence, while assuming charged particle equilibrium from the short-ranged 
delta electrons. At particle-therapy centers air-filled ionization chambers are routinely 
used as a main tool for quality assurance of the delivered beam. Several dosimetry 
protocols for protons have been conceived while the most recent protocol provided by 
the International Atomic Energy Agency (IAEA), TRS-398 |T], sets the standard in 
proton dosimetry today. In addition, TRS-398 also covers dosimetry for ions heavier 
than protons. The protocol uses an absorbed dose-to-water based formalism and relates 
the dose to water D w q to the acquired charge Mq multiplied by a calibration factor 
Nd,w,q and a dimensionless beam quality correction factor /cq,q . The correction factor 
/cq q relates the measured beam quality Q to the beam quality Q used for calibration 
of the dosimeter and it is defined in TRS-398 as 



including the water-to-air stopping-power ratio, Skater/air, the mean energy expended in 
air per ion pair formed, W^, and a perturbation factor Pq/pq which considers effects 
for the specific ionization chamber used. 

As mentioned in TRS-398, calculating the correct beam quality factor in particle 
therapy is complex since it involves knowledge of the entire particle-energy spectrum 
at the point of interest. Instead, TRS-398 proposes a pragmatic approach by 
recommending a fixed value of 1.13 as a generic correction factor for the dosimetry 
of ions heavier than protons based on the analysis by Hartmann et al. [2J, irrespective 
of the particle types and energy spectra which are functions of depth. Accordingly, 
TRS-398 summarizes that the estimated combined standard uncertainty in kg t Q in 
ion beams heavier than protons (about 3%) arises largely from the uncertainty of the 
stopping-power ratio (STPR) (about 2%) and the value for W^ T (about 1.5%) fl This 
has been taken up by Henkner et al. [3] and Geithner et al. [HE] for mono energetic 
carbon beams, and they found out that (i) the STPR is not constant but varies with 
pentration depth, and (ii) it depends strongly on the accuracy of the stopping-power 
data used as input for the calculation. Accordingly, it was concluded that in a clinical 
setting an over or under dosage may occur in the order of a few percent. 

Here, we shall continue the initiated work on STPRs focusing on two objectives. 
First, gaining a sound understanding of the physics determining the STPR, and second, 
exploiting the gained understanding in order to provide results with direct clinical 
relevance which are ready to be applied in clinical practice in a quality assurance setting. 

The deeper insight in the context of STPR is required since STPR strongly depends 
on the stopping-power data which are used for its determination. The problem is, 

| Although the study of the value for is beyond the scope of this work it shall be mentioned that 
the estimated uncertainty for W a i T calls for a detailed investigation. 





(1) 
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however, that the stopping-power data currently recommended by the International 
Commission on Radiation Units and Measurements (ICRU) possess some intrinsic 
inconsistencies. In contrast to an accurate but purely numerical calculation of STPRs, 
a sound understanding of the relevant physics allows for conclusions independent of the 
employed set of stopping-power data. It is also a prerequisit for an analytical description 
of the STPR. It should be emphasized that, in contrast to earlier work on STPR [3H5], 
the present study also considers the recently revised ICRU 73 [5] stopping-power data 
for ions heavier than helium on water. These data replace the ones originally published 
in ICRU 73 which led to a still ongoing discussion on stopping powers for water targets 
(see, e.g., [7]) especially in view of recent measurements by Schardt et al. [8]. 

TRS-398 explicitly states that the STPR for water-to- air, tSwatcr/ain should be 
obtained by averaging over the complete spectra of particles present. And consequently, 
this requirement was considered to be an important limitation in the case of heavy 
charged particles, where the determination of all possible particle spectra was assumed 
to be a considerable undertaking. This was certainly the case a decade ago and may still 
be true from the point of view of dose determination for routinely quality assurance. 
Nowadays, however, the determination of complete spectra of particles can be achieved 
conveniently and with high accuracy by applying Monte Carlo transportation codes 
exploiting the commonly available computer power. These codes are in general valuable 
in predicting radiation fields of ions in tissues and are in particular useful in hadron 
therapy for the simulation of ion transport. The most common codes in particle therapy 
with ions heavier than protons are Geant4 [9], FLUKA [ID], PHITS [H], MCNPX [12], 
and SHIELD-HIT [T3~irH] . all taking into account the atomic interaction of the ions 
with the target medium as well as the nuclear interaction. It is the former interaction 
which mainly determines the energy loss of the incident ions and therefore the stopping 
power, while the latter interaction is responsible for fragmentation and therefore for the 
production of secondary-particle spectra. 

Initial studies on STPRs relevant for dosimetry in radiation therapy with ions 
heavier than protons were performed without Monte Carlo calculations ignoring the 
influence of the secondary particle spectrum (e.g. [151417] as presented in TRS-398 
[1]). Calculations exploiting the capabilities of Monte Carlo codes were performed 
with SHIELD-HIT but exclusively carbon-ion fields [3HS] were studied. However, the 
dependence of the STPR on different ion species is of intrest since a number of facilities 
world-wide (e.g., NIRS and HIT) are equipped with radiation fields which cover a 
broader range of ions than merely protons and carbon ions. Furthermore, it was recently 
argued that ions heavier than carbon may play an important role in the near future 
concerning the radiation therapy of radio-resistent tumors [18] . Consequently, a large 
variety of ion species, namely, H, He, Li, C, N, O, Ne, Si, and Ar are considered here - 
all accessable either for clinical radiation therapy (up till O and Ne at HIT and NIRS, 
respectively) or for in vitro radiobiology experiments. 

Despite their obvious relevance in medical application, so far, STPRs for spread-out 
Bragg peaks (SOBPs) for ions heavier than protons have been discussed only scarcely 
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in literature, namely, by Henkner et al. [3]. In the case of proton beams more detailed 
efforts have been performed, e.g., by Palmans et al. [H?|I2U] and earlier already by Medin 
and Andreo (21]. Henkner et al, who considered carbon ions, outlined in the conclusions 
of [3] that a more detailed analysis of STPRs for SOBPs is clearly needed since their 
statements were only based on the analysis of a single physically optimized SOBP using 
one set of stopping-power data. Consequently, one focus of this work should be a 
systematic study of the STPR for SOBPs, both physically and biologically optimized, 
of different widths and practical ranges leading to an analytic expression for the STPR. 

This paper is organized as follows: First, the physics relevant for the STPR and the 
employed methods are discussed. Furthermore, analytical expressions for the average 
energy of the primary ions and STPRs are proposed. Subsequently, the results for the 
water-to-air STPR for pristine as well as SOBPs are presented and compared to the 
proposed analytical expressions. The following discussion concentrates on three issues, 
namely, the influence of the stopping-power data on the STPR, the dependence of the 
STPR on the ion energy, and STPRs for SOBP. 

2. Materials and Methods 

For all our calculations we used the Monte Carlo particle transport code SHIELD- 
HIT [HUH], based upon the most recent version SHIELD-HIT08 [22]. A number 
of improvements and new functionalities were added to SHIELD-HIT08, documented 
in [23], finally resulting in SHIELD-HIT10A [21]. Here, only the relevant changes are 
reported. First, there is now the possibility of directly scoring the STPR of any media, 
described in detail in section 12.41 Apart from this, raster scan files generated by the 
treatment planning software TRiP [25|l26] can now be read by SHIELD-HIT in order to 
recalculate SOBPs. In this study, we present calculations from four single field carbon 
ion SOBPs, listed in Table [TJ The width of the SOBP is defined as usual by the width 
in which the dose is above 95% percent pQ. All SOBPs are 3-dimensional dose cubes 
with equal side lengths. The resulting raster-scan file describes the needed amount of 
particles for each raster point and for each energy slice providing the necessary input for 
SHIELD-HIT to generate the radiation field for the SOBP. A ripple filter implementation 
based on the design described by Weber et al. [27] is added to SHIELD-HIT in a similar 
way as specified by Bassler et al. in [2S], in order to produce flat SOBPs. 

The practical range, R p , is defined for protons as the depth at which the absorbed 
dose beyond the Bragg peak or SOBP falls to 10% of its maximum value pp. However, 
for ions heavier than protons this definition of R p is not feasible due to the pronounced 
dose tail of secondary particles. Therefore, the depth at which the absorbed dose beyond 
the Bragg peak or SOBP decreases to 50% of its maximum value is proposed and used 
here for ions heavier than hydrogen, i.e., z > 1. Also other definitions of R p have been 
used before as discussed in [22] • The residual range _R rcs at a depth d is than defined as 

-Rres = R p — d (2) 

and the measurement depth <i re f at the middle of the SOBP in accord with TRS-398 p]. 
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2.1. Stopping powers and mean excitation energy 

Stopping power S is defined as the average energy change dE of a particle per unit 
length dl in a medium. At high energies, that is about from 10 MeV/u up to 1 GeV/uJ§| 
the mean energy loss of a charged particle to atomic electrons is well approximated by 
Bethe's original theory [30|l3T] which treats the electromagnetic interaction in first-order 
quantum perturbation theory. At lower energies, however, additional higher-order terms 
are required in order to reproduce experimental results. The transition from the regime 
of quantum perturbation theory to the one permitting a classical treatment is described 
in Bohr's distinguished survey paper [52] . 

A widespread formulation of Bethe's theory summarizing all terms of the lowest- 
order stopping number L was proposed by Fano [33] 

S 47re 4 1ZJ, 2m e v 2 1 n2 C 51 

-p=^uA Z [ ln ^- + ln T^- /3 -Z-2\ ' (3) 

In Eq. (j3J), p is the density of the medium, m e the electron mass, e and u are the elemental 
units of electric charge and atomic mass, respectively, Z and A are the atomic number 
and the relative atomic mass of the target medium, respectively, v and z are the velocity 
and the charge of the projectile, and /? = v/c where c is the velocity of light in vacuum. 
The mean excitation energy of the target medium is denoted by I, while C/Z and 
5/2 are the shell corrections and the density-effect correction, respectively. The second 
and third term in the square brackets containing (3 originate from Bethe's relativistic 
extension [31] and are often referred to as relativistic corrections. The expression in Eq. 
(JSJ) is consistent with the first term L of the stopping number L in ICRU report 49 [T7] . 
For low energies the description of the stopping powers becomes more complicated and 
higher-order terms of the stopping number L have to be taken into account in order to 
correct for a number of different effects, such as the Barkas and the Bloch correction, L\ 
and L/2, respectively. An effective description for the energy regime below the stopping- 
power maximum was provided by Lindhard and Scharff [31] assuming a rise of the 
stopping power which is proportional to the square root of the particle energy. 

§ The energy regime from 10 MeV/u up to 1 GeV/u, corresponds according to the revised tables for 
water in ICRU 73 [S], for carbon ions to a range from 0.0427 cm up to 108.6 cm. 



Table 1. Specifications of four spread out Bragg peaks (SOBPs) for carbon ions in 
water. Given are the width along the beam axis, the practical range i? p , and whether 
the SOBP is optimized for a homogeneous physical dose or relative biological dose. 
The optimization was performed by the treatment planning program TRiP (251126) . 

SOBP Width (mm) R p (mm) Optimization 

a ~~ 50 220 physical 

b 80 168 physical 

c 50 150 physical 

d 100 153 biological 
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Table 2. Specifications for 6 sets of stopping-power data used in this work. The 
stopping-power data for the first three sets are determined internally by SHIELD-HIT 
(cf. Sec. 12. 2[) using the given values for I W ater and I a i r . while those for the sets 4 to 6 
are directly read by SHIELD-HIT as text files in tabulated form. For the latter, two 
different tables per set are used distinguishing between the lightest (H and He) and 
heavier ions. The table specifies for each set its number, / wa ter an d -fair in eV, the 
range of ions for which these data are applied, references, and if adequate additional 
comments. Further explanations can be found in the text. 



SHIELD-HIT calculates stopping-power data using /-values 
Set # /water -fair ion range Reference Comments 



1 


78 


82.8 


z > 1 


ICRU 73 


using revised / watcr [32] 


2 


75 


85.7 


z > 1 


ICRU 49 pi] 




3 


80.8 


85.7 


z > 1 


Henkner et al. [3] 






SHIELD-HIT directly 


uses tabulated stopping-power data 


Set # 


/water 


/air 


ion range 


Reference 


Comments 


4 


78 


82.8 


z > 2 


ICRU 73 [61136] 


revised data for water [36] 




75 


85.7 


z < 2 


ICRU 49 [13 




5 


67.2 


82.8 


z > 2 


ICRU 73 [6] 


only original data for water 




75 


85.7 


z < 2 


ICRU 49 pi] 




6 


75 


85.7 


z = 1 


ICRU 49 PH 






75 


85.7 


z > 1 


MS TAR [37] 


charge scaling of ICRU 49 



The mean excitation energy, /, is a property of the medium which enters logarithmi- 
cally in the stopping formula Eq. fl3]), and is responsible for most of the target material 
dependence of the stopping-power. It is, on the other hand, completely independent 
of the properties of the projectile. According to Eq. fl3]), a larger /-value results in a 
smaller stopping power and consequently in a larger range of an ion in the medium. 
The /-values in the ICRU report 49 [17] for protons and alpha particles (retained from 
ICRU report 37 for electrons [35]) were mainly taken from measurements. In ICRU 
report 73 [6], however, the /-values are mostly determined theoretically As a result 
different /-values for the same material are recommended in ICRU reports 49 and 73. 
Obviously, this is inconsistent, since an /-value should not depend on the projectile. 
The differences existing between the ICRU reports highlight that the accuracy of the 
current employed methods to determine stopping-power data have still to be improved 
in order to provide a consistent target description. 

2. 2. Stopping powers in SHIELD-HIT 

In the current implementation of SHIELD-HIT the compilation of required stopping- 
power data can be done in two ways which can be chosen independently for each target 
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medium. First, stopping-power data can be calculated internally by SHIELD-HIT using 
a modified Bethe formula at high energies and a Lindhard-Scharff description [34] at low 
energies for any kind of material composition using the corresponding material-specific 
values for I, Z, and A as discussed before in Sec. 12.11 Second, an arbitrary stopping- 
power table may be read in as a formatted text file allowing for the use of, in principle, 
any stopping-power data which can be provided in electronic form. In this work the 
common open source library UbdEdx |38, 39] which is available online is applied in order 
to provide tabulated data in formatted form from the ICRU reports 49 [TTJHU] and 
73 PES] as well as MS TAR [37]. 

The Bethe formula used by SHIELD-HIT is similar to the formulation in Eq. ([3]). 
But, so far no shell corrections C/Z have been considered. These are known to be 
most relevant for low energies where, however, the Lindhard-Scharff description is used 
instead in SHIELD-HIT. Furthermore, it was demonstrated that for low energies (about 
1 MeV/u) the accuracy of stopping-power data is insignificant for particle therapy [4T] . 
The same argument holds for the higher-order term L\. Additionally, the Bethe formula 
is modified in order to allow for electron capture (significant for low energies) by using 
an effective energy-dependent scaling of the projectile charge z by Hubert et al. |42j . 
Currently, relativistic corrections proposed by Lindhard and S0rensen [43] are still 
missing in SHIELD-HIT. Their importance increases for heavy ions with large nuclei 
which cannot be approximated as point-like particles. Although their relevance for 
particle therapy should be studied no significant impact has been expected so far. 

Due to existing inconsistencies in the stopping-power data recommended by ICRU 
- discussed in Sec. 12.11 — different sets of stopping powers are used in this work, all 
listed in Table [2j Thereby, sets 1 and 2 as well as sets 4 and 5 are directly related to ICRU 
reports. For comparison, in set 3 the preferred /-values of Henkner et al. [3] are used 
while set 6 employes the frequently used data provided by MSTAR [37J. The intended 
purpose of the sets 1 and 2 is the attempt to describe the target media consistently 
with only one /-value for all ions, both with z < 2 as well as z > 2, applying SHIELD- 
HIT's internal routine to determine the stopping power. Accordingly, set 1 uses only 
the /-values from ICRU report 73, / 73 , (the revised value for water, / wate r = 78 eV, 
was very recently published in the erratum to ICRU 73 [36] ) while only /-values from 
ICRU 49, Z49 are used in set 2. The motivation for sets 4 and 5, on the other hand, is 
the direct application of the recommended tabulated data which can be found in ICRU 
reports 49 and 73 for ions with z < 2 and z > 2, respectively. While set 4 uses the 
recently revised stopping-power data for heavy-ions on water, set 5 uses, for comparison 
to earlier studies of the STPR, the water data as originally published in ICRU 73. Note, 
the recently revised data from ICRU 73 [36] were not employed by Henkner et al. [3]. 
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2.3. Stopping-power ratio 

The stopping-power ratio S a /b between medium a and medium b is (cf. TRS-39s[j] pQ) 
given as a particle fluence weighted average over all primary and secondary particles. It 
is determined by calculating the dose ratio via track- length fluence <§> a ^(E) of particle % 
in medium a as function of particle energy E and mass stopping power Si(E)/p 



In Eq. (j2J) numerator and denominator are equal except for that the mass stopping power 
of medium a enters in the numerator and of medium b in the denominator. An energy 
cutoff E min > may originate, e.g., from the chamber geometry. The contribution of 
"track-ends" to the total dose deposition and to the corresponding STPR was studied 
in [H| . There it was concluded that they are not of relevance for light-ion dosimetry 
which is in contrast to electrons, where the contribution to the total deposited dose can 
be between 6% and 8% [55] . 

In contrast to the correct definition for the STPR of an ion field in Eq. (J4]), the 
ratio of stopping powers for media a and b for one particle species of energy E, 



has often been considered as an approximation to the STPR, e.g. in [UEIE]- The right 
hand side of Eq. (J5J) is expressed by Bethe's stopping formula as given in Eq. (J3J) but 
omitting corrections. Note, the ratio in Eq. fl5]), which considers only one particle species, 
is a function of the particle energy E in contrast to the STPR in Eq. (J3J) which has a 
spatial dependence and takes the full energy spectra of all particles into account. 

2.1 Scoring of STPR in SHIELD-HIT 

STPRs have already been obtained with SHIELD-HIT before and only the 

conceptual improvements in this work are discussed in the following. The concept of 
virtual scoring has been introduced which now allows for a parallel detector geometry 
independent of any physical geometry. Therefore, there is no longer a need for 
introducing artificial physical geometries which lead to additional region boundaries. 
Furthermore, the STPRs are now determined on-line, that is, during the transport of 
the particles. An on-line calculation has the advantage that possible influence on the 
result due to the number and size of the energy scoring and energy spacing is avoided. 
Additionally, higher accuracy in scoring of tracks-ends can be achieved in principle. 

The detector for the STPR resembles Eq. (j2J) and is implemented in the following 
way. When a particle traverses a bin of the STPR detector its track-length fluence 
within the bin is scored and directly multiplied with (S/p) a of the medium a in which 

|| In IAEA TRS-398 only the water-to-air STPR is explicitly denned, i.e., a=water and 6=air. However, 
this definition is also useful for other media combinations. 




Ejgj^m (Sj(E)/p) a dE 
ZJZ^aAE) (Si{E!)/p) b dE- 



(4) 



(S(E)/p) a _ (i) a \n[2m eV yi a } 
(S(E)/p) b (f > fe H2m e vyi b ] ' 



(5) 
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the particle moves for the energy (E in + E out )/2. E in and E out are the energies of 
the particle when it enters and leaves the bin, respectively. Additionally, the same 
track-length ffuence is multiplied with (S/p)b of the same particle. Both quantities are 
summed up individually including all particles passing the bin. After a full Monte Carlo 
transport simulation the two sums are divided yielding the STPR for this bin. 

In this work a transport cutoff of 0.025 MeV/u is used by SHIELD-HIT which means 
that all particle tracks end once the particle energy becomes smaller^ Consequently, 
the lower limit for the integration in Eq. (j4j) is given by -Emm = 0.025 MeV/u having an 
influence on the STPR of less than 0.00015% [3]. A recent review article [15] discusses 
in some detail the impact of electrons in fast ion-atom collisions with respect to hadron 
therapy as well as the possibility to extend SHIELD-HIT in a way that also electron 
tracks are considered. This would allow for studies of the microscopic energy distribution 
in the target medium. Tracking of delta-electrons has for example been performed with 
Geant4 |3B] and a study comparing to the present work might be of interest. 

2. 5. Analytic expression for the average ion energy 

The stopping-power formula as presented in Eq. for a specified combination of 
projectile and target is primarily a function of the projectile's kinetic energy which 
decreases during the passage through the target medium due to the energy loss. In 
order to determine the average energy of the projectiles as a function of depth a full 
simulation of the particle transport has to be performed. This comprises the slowing 
down caused by all relevant energy-loss mechanisms including elastic as well as nonelastic 
interactions [22]. Thereby, nonelastic nuclear reactions produce a spectrum of particles 
with each particle having an individual energy distribution which is furthermore a 
function of the position in the medium. Consequently, it would be highly desirable to 
have a simple, though approximate, analytical expression E(d) for the average energy 
of the primary particles with initial energy Eq as function of penetration depth d. 

Starting with the Bethe formula, but assuming first that the expression in the square 
brackets of Eq. (J3J is independent of energy, E(d) can easily be expressed analytically, 



where R p is the practical range. R p depends in general on the ion species, E , and the 
target material. For energies relevant in particle therapy R p can often be approximated 
by Rq obtained with the continuous slowing down approximation (CSDA). 

In order to account for the correct energy dependence of the Bethe formula as well 
as nonelastic collisions one has to allow for a more general power-law relation, 



This is consistent with ICRU 73 [6] where the range tables for liquid water show the average path 
length travelled for slowing-down from initial energy E to Eq = 0.025 MeV/u. 




(6) 




(7) 
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with an exponent k. Different values for k are suggested in the literature while Kempe 
and Brahme [2H] proposed the use of a dimensionless transport parameter k = E q /R q Sq 
with So = S(Eq). A simple value of about k = 1.7 fits the calculations performed with 
SHIELD-HIT being also compatible with and is therefore used in this study. 



2. 6. Analytic expression for STPR 

In order to derive an analytic, though approximate, expression of the STPR as a function 
of the depth d for two media a and b, the approximation to the average energy in Eq. (|7j) 
can be used together with the ratio of stopping powers given in Eq. (j^J). Utilizing the 
non-relativistic relation v 2 = 2E/m p between the particle velocity v and and its kinetic 
energy E, where m p is the proton mass, one obtains the expression 

§ (A (l)q HEo/Q + c{d) 

£(a/b) (d) - jzV T-rErTFT—7yJ\ ^ 



, A 



) b \n[E /I h ] + C(d) 



where 



C(d) = - In 
k 



d 

1 ~r 



v J 



6.1291 (9) 



and In [4 m e /m p ] = —6.1291 have been used. Similar as in Eq. fll]) the numerator and 
denominator in Eq. (jSJ) equal except for the different /-values and (Z/A) ratios. It 
should be mentioned that in order to keep the expression for S^a/b) as simple as possible 
its derivation has been performed without relativistic kinematics which are in principal 
of relevance for the highest energies used in particle therapy. Finally, the expression in 
Eq. (jHJ) should explicitly be formulated for the water-to-air STPR 

5 {waWair )(rf) - 1.1H95 HEo/hir] +1/Jfc i n[1 _ d/jRp] + 7. 686 3 ' ^ 10 ) 

being the most relevant case for dosimetry in particle therapy, with E and / in units of 
MeV/u and eV, respectively. For convenience, constants are expressed in numbers, i.e., 
ln[4 x 10 6 m e /m p ] = 7.6863 and (f ) / <f ) = 0.555076/0.499189 = 1.11195 [35]. 



3. Results 



In what follows, results for STPRs of pristine Bragg peaks and SOBPs are presented. 
Furthermore, analytic expressions are compared to the numerical results obtained with 
SHIELD-HIT. For simplicity, only results for ion beams in water are considered, and 
STPRs are presented exclusively for water-to-air. Furthermore, only stopping-power 
set 1 of Table [2] is used in order to ensure consistent results except for Fig. QJb) which 
demonstrates the dependence of Skater/air on the choice of the stopping-power set. 



3.1. Pristine Bragg peaks 

A comparison between the calculated STPR for water-to-air of 270 MeV/u carbon ions 
as a function of depth in water and the corresponding depth-dose distribution is shown 
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110 120 130 140 R p 150 20 40 60 80 100 120 140 160 180 

Depth (mm) Depth (mm) 



Figure 1. 270 MeV/u carbon ion beam in water, (a) Depth-dose distribution 
(line) and stopping-power ratio Skater/air (symbol) as a function of depth around the 
Bragg peak. The stopping-power set 1 is used, (b) S watcr / air as function of depth 
obtained with the six stopping-power sets 1 - 6 as specified in Table [5] and the value 
1.13 recommended by IAEA in TRS-398 [TJ. 

in Fig. QJa) focusing on the vicinity of the Bragg peak. The maximum of the STPR 
almost coincides with the practical range R p = 144.6 mm and therefore appears to be 
at a larger depth than the dose maximum. The width of the STPR peak is considerably 
smaller than that of the dose curve. The determined height of the STPR depends to 
some extend on the finite spatial resolution along the beam axis. A coarse resolution 
leads to spatial averaging and accordingly to a lower peak height of the STPR. Note, 
the position of the Bragg peak and therefore R p is only influenced by the choice of 
stopping-power data for water. 

The influence on the STPR due to the use of different stopping-power data is 
demonstrated in Fig. QJb) for a 270 MeV/u carbon pencil beam by using all stopping- 
power sets 1 to 6 of Table |2j In the plateau region, the deviations of the sets 1 to 6 
from the value 1.13 recommended by the IAEA in the TRS-398 |Ij is within 1%. Set 
2, which employs I49 for water and air, differs the least. This is consistent since the 
recommendation in TRS-398 is based on the stopping-power data provided by ICRU 
report 49. In contrast to the recommended value 1.13, none of the calculated STPR 
curves is constant. However, the relative increase in the plateau region up to a depth 
of 130 mm (-R res ~ 15mm) is moderate and of the order of approximately 0.2% to 
0.3%. For all sets of stopping powers, except for set 4, an increase of the STPR can 
be observed in the vicinity of the Bragg peak. Set 4 on the other hand shows a dip. 
This dip originates mostly from the carbon ions and can therefore be attributed to the 
tabulated data provided by in the revision of ICRU 73. 

The two STPR curves calculated with set 1 and set 3 — the latter was used in [3] 
- lie virtually on top of each other in Fig. QJb) although the /-values of set 1 and set 
3 differ notably. This can be explained by the differences of the /-values for water and 
air which are very similar with 4.8 eV and 4.9 eV for set 1 and set 3, respectively. As 
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Figure 2. Stopping-power ratio for water-to- air S^to-c / air a s a function of i? rcs in 
water. The stopping-power set 1 specified in Table [2] is used, (a) Carbon ions with 
different initial energies ranging from 70 MeV/u to 400 MeV/u are compared, (b) 
Beams of different ions relevant for particle therapy are compared: H, He, Li, C, N, O, 
Ne, Si, and Ar. Ion beam energies ranging from 130 MeV/u to 400 MeV/u have been 
chosen in order to achieve comparable penetration depths. 



expected, the S water / air obtained with the tabulated data from ICRU 73, that is set 4, 
agrees with the STPR curve obtained with / 73 , set 1, in the plateau region. On the other 
hand, the STPR curves for the sets 1 and 4 deviate around and beyond R p . The use 
of the stopping-power set 5 by Henkner et al. in j3] resulted in an unphysical minimum 
of the STPR in the plateau region. In the present study, however, no minimum of the 
^water/air curve obtained with set 5 can be observed. 

The influence of different initial energies on the STPR for carbon ion beams as 
function of the residual range R rcs is shown in Fig. [2(a). The STPR curves as function 
of R res are almost identical in the plateau region and therefore independent of the initial 
energy. Around and beyond the Bragg peak the curves are still alike, though lower 
initial energies lead in general to higher STPRs values. 

Stoppin g-pOWer ratio ^water/air SiS a function of R res for different ion beams relevant 
for particle therapy: H, He, Li, C, N, O, Ne, Si, and Ar are presented in Fig. [2(b). 
Different beam energies for the individual ions have been chosen, ranging from 130 
MeV/u to 400 MeV/u, in order to achieve comparable penetration depths. In the 
plateau region the STPRs for the different ion beams all share the same qualitative 
behavior which has already been discussed before in the context of carbon ions, cf. Figs. 
H(b) and [5(a). A comparison among the ions yields that decreasing STPRs at a given 
R res occur for increasing atomic numbers z. The decrease becomes less pronounced for 
larger z. An exception from this trend is found for H and He ions with very similar 
<Swater/air m the plateau. The relative difference of STPRs between the lightest and 
heaviest ion z = 1 and z = 18, respectively, is rather constantly about 0.15%. 

Around and beyond R p (R Tes ^ 0), the STPR seems to be larger for ions with larger 
z. Although, some dependence also might originate from the different initial energies 
of the ion beams, as has been observed in Fig. [2(a). While the STPR for H ions is 



Analytical expressions for water-to-air stopping-power ratios 



13 




1 1 1 1 1 1 1 1 




TRS-398-^ 


(b) 






$ 




5 Oram, phy ; 






80mm, phy _ 


1 




50mm, phy - 






100mm, bio - 








FIT 

" 









*" 



20 40 60 80 100 120 140 160 180 200 220 240 

Depth (mm) 



■20 o 



20 40 60 80 100 120 140 160 
Residual range (mm) 



Figure 3. Different SOBPs for carbon ions in water as specified in Table [T] The 
stopping-power set 1 specified in Table [2] is used, (a) Depth-dose distributions (line) 



and STPRs S. 



watcr/a 



(symbol) as a function of depth for the biologically optimized 



SOBP d. (b) STPRs S, 



watcr/a 



as a function of the residual range R Tes for all four 



SOBPs a, b, c, and d. Also shown is a simple fit for set 1 proposed in Eq. (|TT|) . 



nearly identical to that for He ions in the plateau region, differences to all other ions 
with z > 1 are obvious for i? res < 0. However, as discussed in Sec. 121 the definition 
of the practical range R p for protons differs from that used for ions with z > 1 which 
influences — according to Eq. (j2J) — also the residual range i? res . 

3.2. Spread-out Bragg peaks 

Four different SOBPs obtained with carbon ion fields are considered with the stopping- 
power set 1 in this study and their properties are listed in Table [TJ Figure [3(a) 
displays the depth-dose distribution (line) and the corresponding STPR (symbol) of 
the biologically optimized SOBP d. The proximal start of the SOBP region can be 
recognized in STPR curve as a significant increase. Towards the distal end of the 
SOBP the STPR reveals an exponential increase. As for pristine peaks in Fig. [1(a) 
the maximum of the STPR curve is reached close to R p . A sharp fall off of the STPR 
occurs for depths d > R p which finally results in a constant STPR. Note, this qualitative 
description is valid for all SOBPs a to d and is therefore independent of the specific form 
of the SOBP or whether it is physically or biologically optimized. 

The STPRs of all four SOBPs a, b, c, and d are compared in Fig. [3Jb) as a function 
of Rres- The SOBPs a and c share the same width but differ in R p . Therefore, their 
STPRs as a function of R rcs are nearly identical. The STPR curves for b and d show a 
very similar behavior as a and c but are extended to 80 mm and 100 mm, respectively, 
according to the larger width of their SOBP region. 

It should be mentioned that for SOBP b a number of ripples in the STPR can 
be observed in the proximal SOBP region. They originate from the optimization 
program TRiP which assumes treatment at the SIS accelerator at Gesellschaft fiir 
Schwerionenforschung (GSI), Darmstadt, Germany. The SIS provides finite energy steps 
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Figure 4. SOBP c for carbon ions in water as specified in Table [TJ The stopping- 
power set f specified in Table [2] is used. Dose distribution (line) and STPR <!5water/air 
(symbol) as a function of transverse position at a depth d IC { = 150 mm. 
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Figure 5. Stopping-power ratio iS wa ter/aii f° r a 270 MeV/u carbon ion beam in 
water. The stopping-power set 1 specified in Table [5] is used. The STPR obtained 
with the full particle spectrum is compared to two different approximations: STPR 
determined (i) by considering only the contribution from carbon ions, and (ii) with 
the analytic expression proposed in Eq. (|10[) . (a) Absolute Swatcr/air as a function of 
depth, (b) Relative difference between the STPR for the full particle spectrum and 
the two approximations to the STPRs presented in (a) as a function of residual range. 

which are more coarse at the lowest energy part when no bolus is applied. 

The calculated STPR and dose transverse to the beam axis at the reference depth 
d TC f = 150 mm (defined as middle of SOBP [T]) is displayed in Fig. H]for the SOBP c. 
The STPR is perfectly constant within the full extension of the SOBP transverse to the 
beam axis. A moderate increase of the STPR occurs outside the SOBP. 

3.3. Analytical description of STPR 

The purpose of this section is to relate an analytical description of the STPR to the 
numerically obtained STPR. Thereby, it is important to keep in mind that according to 
the results obtained so far the STPR (a) depends primarily on R Tes , that is average ion 
energy, and (b) is qualitatively independent of the ion species. 
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The STPR Skater/air for a 270 MeV/u carbon ion beam using the stopping-power set 
1, as shown before in Fig. [T^b), which considers the full particle spectrum is compared 
in Fig. [5] to two approximations of the STPR: (i) STPR obtained with SHIELD-HIT 
but ignoring the influence of produced fragments on Skater/air by only considering the 
STPR resulting from carbon ions and (ii) the analytic expression in Eq. ffTOl which 
approximates the STPR with a ratio of stopping powers of the primary ions with an 
average energy depending on i? rcs . Figure E^a) clearly shows that especially in the 
plateau region the absolute difference between these three curves is small. Therefore, 
Fig. 0(b) additionally shows the relative difference \S — S appr | / S fu11 between the 
STPR for the full particle spectrum, determined according to Eq. 03]), and the two 
approximations of the STPR. Both approximations reproduce S iul1 within 0.02% in the 
whole plateau region. Around R p the difference can be as large as 1% and they cannot 
be applied beyond R p since both approximations are based on the primary particles. 

In principle, the analytical expression for pristine peaks can also be of use for SOBP. 
However, since the exact weights of the superposed peaks are not always known a simple 
fit as a function of d may be proposed 

Swatcr/air^) = a + ft exp[7(i? p - 6)\ + 6 (Rp ~ d) (11) 

which approximates the STPR within the SOBP region and is also shown in Fig. [3]^b). 
The values for the four parameters in Eq. (II ip are 

a = 1.12205, (3 = 4.0044E-03, 7 = -0.241, and 5 = -2.0238E-05 . 
These parameters depend on the stopping-power set and slightly on the ion species [39J . 
Outside the SOBP region the STPR might be approximated in a similar way as it is 
done for pristine peaks considering the average energy of the primary ions at d which 
is, however, different from the case of pristine peaks. From Fig. &b) it can be observed 
that the fit function proposed in Eq. ffTTl) clearly is in acceptable agreement with the 
Swatcr/air curves of the four different SOBPs. 

4. Discussion 

4-1. Influence of the inconsistency of ICRU stopping-power data on the STPR 

It is obvious that the STPR strongly depends on the stopping-power data which are 
used as an input for its determination. This applies for R p and accordingly the position 
of the STPR maximum. Therefore, the present work focuses mainly on stopping-power 
data for water and air which are recommended by ICRU reports. However, even if one 
tries to follow these recommendations certain inconsistencies still remain and different 
sets of stopping-power data — as listed in Table [2] — can be deduced. In the case that 
the tabulated data provided by ICRU are used directly, set 4, physical quantities of 
the target media, e.g., the /-value are depended on whether the medium interacts with 
protons and helium ions or with heavier ions. Although not explicitly recommended, 
another approach, which is more consistent from a physical point of view, is to use the 
/-values of only one of the ICRU reports together with an appropriate stopping formula 
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for all ions, i.e., with z < 2 as well as z > 2. This is done here for set 1 and set 2 using 
I 73 and J 4 9, respectivelyE] 

First, the STPRs in the plateau region (d < R p ) are discussed. Figure [TJb) reveals 
that the STPR Skater/air using J 4g for all ions, set 2, agrees best with a constant value 
of 1.13 as recommended in TRS-398 which is plausible since TRS-398 is based on the 
data of ICRU 49. Skater/air obtained with 7 73 , set 1, is about 1% smaller compared to 
that obtained with I^ 9 , set 2 but agrees nicely with the tabulated data recommended 
by ICRU, set 4. Second, around R p the STPR curves of all data sets show a distinct 
maximum except for that of set 4 for which a minimum can be observed. This minimum 
does not originate from differences in the target descriptions due to the use of ICRU 
49 and 73 in set 4 but is exclusively caused by the low-energy ratio of stopping powers 
taken from ICRU 73. Third, beyond the peak (d > R p ) S wate r/air is rather constant for 
all stopping-power sets with a consistent description of the target media. For set 4 and 
set 5 the faster decline of the heavier ions relative to the lighter ones beyond R p leading 
to a transition from ICRU 73 to 49 is clearly revealed byQJb). Therein, the latter two 
curves finally converge to S water /air obtained with / 49 , set 2. 

The use of the out-dated standard, that is stopping tables from ICRU 49 and 
73 without the revisions of ICRU 73 for water, set 5, is not advised. Set 5 yields a 
STPR which is in the plateau region about 1% and 2% larger than the value 1.13 
recommended by TRS-398 and the values obtained with revisions of ICRU 73 for water, 
set 4, respectively. It should be mentioned, that the results for set 5 in Fig. [TJ(b) do not 
show any unphysical minimum in the plateau region as was observed before in [3JH]. 
There, the unphysical structure was attributed to the use of different stopping-power 
data from ICRU 49 and 73 for different ions. However, in the present work it was 
possible to reproduce the exact shape of the curve as shown in Fig. 2 of [3] by using 
only three digits of the stopping-power data of ICRU 49 and 73 instead of the four 
digits provided by the tables and therefore to contradict the earlier statement. We can 
conclude that the reason for the unexpected behavior observed in the corresponding 
Fig. 2 of [3] is simply due to errors resulting from a too coarse rounding in the applied 
stopping-power list and are not caused by a combined use of ICRU 49 as well as 73. 



4-2. Dependence of the STPR on the average ion energy 

One quintessence of this work is that the STPR for a given set of stopping-power data is 
mostly determined by the average energy of the ions, rather than their initial energy or 
their charge which is nicely confirmed by Fig. |2] Figure [5] shows furthermore that for d < 
R p the STPR is completely dominated by the STPR of the primary ion species and the 

+ A further option is to correct for the /-value in one of the two recommended sets of ICRU tables 
leading to a more consistent description of the target media. According to Eqs. (|3|) this could be done in 
first order using the term 0.307075 (z 2 Z)/{0 2 A) \n\I49 / for correcting ICRU report 49 or its negative 
value for correcting ICRU report 73 [48]. This approach is of course not applicable in the low-energy 
regime where, on the other hand, the ICRU tables anyhow provide a limited accuracy only. However, 
this option has not been pursued in this work. 
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relative deviation is of the order of 0.02%. In a next step, the STPR of the primary ions 
can be nearly exactly reproduced using the ratio of stopping powers, as expressed in Eq. 
(EJ), for the average energy of the primary ions at a depth d. Since the average energy 
can be rather accurately expressed as a function of d, in a final step this energy function 
is used together with Eq. fl5]) to formulate the analytical expression of the STPR in Eq. 
(II Op . The analytical expression reproduces the STPR of the primary ions very well and 
deviates accordingly also in the order of 0.02% from the correct STPR obtained with 
the complete particle spectrum. 

Note, the key advantage of the proposed analytical expression is its flexibility. It is 
not restricted to a specific set of stopping-power data. Consequently, it could be easily 
adopted to any new recommendation by ICRU. It is not restricted to specific primary 
ions and only their average energy as a function of depth is required. It is not restricted 
to a specific combination of target materials such as water and air being in the focus 
of the present work. For example, it can be straightaway used for STPR for water to 
tissue and air to tissue which would be of interest when comparing dose to medium with 
dose to water as recently discussed by Paganetti [3H|- 

These findings lead to two central insights. First, in contrast to the presumption 
in TRS-398, the knowledge of the whole particle spectra is not of practical importance 
for the plateau and peak peak regionjf] The STPR can be simply approximated with 
a relative deviation much smaller than the uncertainties of all available stopping-power 
data. On the other hand, the secondary particles are of central importance beyond R v 
where the primary ions are ceased. Second, it is not necessary to study STPRs for all 
ions independently since they can be approximated in the same way as proposed in Eq. 
(JBJ). The small quantitative differences among the STPRs of the various ions shown in 
Fig. Mb) may be explained best with the different average energies at a given R ves . 

4.3. STPRs for SOBPs 

Clinical applications require a relatively uniform dose to be delivered to the volume to 
be treated and for this purpose the ion beam has to be spread out both laterally and in 
depth. With respect to the STPR of SOBPs the gist of this work is that the qualitative 
behavior of the STPR is hardly dependent on the specific spatial form of the SOBP, its 
depth in the medium, and whether it is optimized for a homogeneous physical dose or 
relative biological dose. 

This statement is nicely verified in Fig. [3^b) where Swater/air is displayed as a function 
of i? res for four SOBPs specified in Table [TJ In order to understand the observed uniform 
behavior of the STPR one has to keep in mind that a SOBP is a superposition of Bragg 
peaks of different intensities and R p which is usually constructed from the distal end, i.e. 
R res = 0, toward the proximal start. Consequently, the properties of the SOBP up till a 
residual range R Tes are only weakly influenced by the properties of the SOBP for larger 

* The relevance of the whole particle spectra for determining the STPR increases if for some secondary 
particles different stopping-power data are employed, e.g. ICRU 73 and 49. 
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R res . It is therefore possible and practical to propose a fit to the STPR .Skater/air for the 
SOBP region, as it is done in Eq. (fTTj) for carbon ions and set 1. A general drawback 
of the fit function compared to an analytical expression is that the parameters of the 
former explicitly depend on the ion species and stopping-power data. In analogy to the 
findings for pristine peaks quantitative differences can be expected for ion species other 
than carbon ions. However, a detailed study of a number of other ions is beyond the 
scope of this article and might be addressed elsewhere [3H] ■ The qualitative dependence 
of Skater/ah- on the stopping-power set is similar as discussed before for the pristine peaks. 
Accordingly, the quantitative difference between Skater/air obtained with set 1 and set 4 
is due to the two features observed for set 4, namely, the minimum at R p and the strong 
increase of the STPR beyond R p caused by the ICRU 49 tables for protons and helium. 

In TRS-398 reference conditions for the determination of absorbed dose for ion 
beams are specified. The reference depth d rc { for calibration should be taken at the 
middle of the SOBP, at the center of the target volume. It can be seen in Fig. [5] that 
a positioning error of a dosimeter transverse to the beam axis has no relevance on the 
STPR as long as the position is within the SOBP. This is plausible since the average 
energy of the ions should be same at the same depth. A misalignment along the beam 
axis, on the other hand, may have an influence as seen in Fig. [3j The influence is largest 
for a SOBP with small width, for which the gradient of STPR is largest, and becomes 
smaller for large widths. Therefore, an extended SOBP might be recommended for 
accurate dose measurements in a practical quality- assurance setting. The total variation 
of the STPR along the beam axis observed in Fig. [3]^b) is of the order of 0.8%. 

5. Conclusions 

Calculations of the water-to-air stopping-power ratio (STPR) Skater/air using the Monte 
Carlo transport code SHIELD-HIT10A are performed for different ions in a range of 
1 < z < 18. The STPR is determined on-line considering the track-length fluence 
spectra of all primary and secondary particles as recommended by IAEA in TRS-398. 
In addition to providing accurate quantitative results the focus of this work is put on a 
thorough qualitative understanding of the dependencies of the STPR and the relevance 
for particle therapy. 

STPRs obtained with different sets of stopping-power data recommended by ICRU 
[Hid?], including the very recently revised data for water [35], are compared with the 
value 1.13 recommended for Skater/air in TRS-398 [1] resulting in deviations of the order 
of 1% in the plateau region. The change of the STPR due to the contribution of 
secondary particles is only of the order of 0.02% for pristine peaks in the plateau region 
and up to the Bragg peak. It can be shown that for a given set of stopping-power data 
the STPR at a residual range R res is mostly determined by the average energy of the 
primary ions, rather than their initial energy or their charge z. A convenient analytical 
expression for the STPR as a function of depth in water is proposed for the plateau 
region up to the Bragg peak which deviates in this region by about 0.02% from the 
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obtained results for Skater/air- The most valuable property of the analytical formula 
is its flexibility. It is in principle not restricted to any specific ion, stopping-power 
data, combinations of target media, or initial ion energies. For the case of spread-out 
Bragg peaks (SOBPs) it can be concluded that the qualitative behavior of the STPR 
is hardly dependent on the specific spatial form of the SOBP, its depth in the medium, 
and whether it provides a homogeneous physical dose or relative biological dose. A fit 
function is provided to approximate the STPR within the SOBP region for carbon ions. 

Finally, it can be stated that no further theoretical studies of STPRs heading only 
for higher accuracy are expedient, as long as no consistent set of relevant stopping-power 
data for all ions is recommended, preferably with smaller uncertainties. 
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